Granular segregation under vertical tapping 
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We present extensive Molecular Dynamics simulations on species segregation in a granular mixture 
subject to vertical taps. We discuss how grain properties, e.g., size, density, friction, as well as, 
shaking properties, e.g., amplitude and frequency, affect such a phenomenon. Both Brazil Nut 
Effect (larger particles on the top, BN) and the Reverse Brazil Nut Effect (larger particles on the 
bottom, RBN) are found and we derive the system comprehensive "segregation diagram" and the 
BN to RBN crossover line. We also discuss the role of friction and show that particles which differ 
only for their frictional properties segregate in states depending on the tapping acceleration and 
frequency. 



Granular materials are systems of many particles inter- 
acting via short ranged repulsive and dissipative forces, 
both normal and tangential to the surface of contact. 
They are characterized by an energy scale mgd (of a grain 
of mass to and linear size d in the gravitational filed <?) 
which is many orders of magnitue larger than the thermal 
energy ksT, and are thus named "non thermal" systems. 
These characteristics make difficult the understanding of 
the large variety of counterintuitive phenomena granu- 
lar materials exhibit, which are of great interest both for 
their industrial relevance, and for the theoretical chal- 
lenges posed to physicist and engineers. 

Particularly the phenomenon of size segregation un- 
der vertical vibrations which we consider here, has 
emerged as a real conundrum. Contrary to intuition, 
an originally disordered mixture when subject to ver- 
tical vibrations tends to order: large particles typi- 
cally rise to the top, as smallparticles percolate into 
their voids during shaking Q, 0, 0, H| or move to the 
bottom due to convection mechanisms H, IS 0j giv- 
ing rise to the so called "Brazil Nut Effect" (BN). Dif- 
ferences in particle density also affect size separation 
(see ref.s in [$]) and reverse-BN (RBN), with small 
grains above, can be observed too The picture 

where grain sizes and weights are the parameters ex- 
plaining segre gati on is found, however to be too simple 
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scenario is still missing. 

In correspondence with some existing experiments 
[2ll I22I. Eij ] , here we consider segregation phenomena in 
Molecular Dynamics simulations of tap dynamics: grains 
confined in a box are shaken and after each shake fully 
dissipate their kinetic energy before being shaken again. 
A picture of our model system is given in the left panel 
of Fig. ^ showing the final BN configuration reached by 
an initially disordered mixture shaken with an amplitude 
r = Au> 2 /g = 1 (where oj is the shake frequency, A its 
amplitude and g gravity acceleration, see below). An 
example of the role of the external drive on segregation 
can be appreciated by comparison with the right panel 




FIG. 1: We show a mixture of Ni = 240 large particles of 
diameter Di — 1 cm and density pi — 1.9 g cm -3 (dark gray 
particles) and N s — 360 small particles of diameter D s — 0.8 
cm and density p s = 1.27 g cm" 3 (white particles). It is con- 
tained in a box (whose base is made of other immobile grains, 
light gray particles, see text) and is subject to vertical taps 
with normalized amplitude F. The pictures show two config- 
urations at rest attained at stationarity: interestingly, when 
shaken with T = 1 (left) the system goes into a BN configura- 
tion and when F = 3 (right) it goes in a RBN configuration. 



of Fig. [^showing the final RBN configuration reached by 
the same mixture when shaken at T = 3. 

We show below how grain properties, e.g., size, density, 
friction, as well as the external forcing, e.g., shaking am- 
plitude and frequency, affect the process and derive for 
the first time a comprehensive non trivial "segregation 
diagram" . The richness of such a diagram is not cap- 
tured by current theoretical approaches 0, E, 0] and 
calls for new theoretical and experimental investigations. 

Simulations - We make soft-core molecular dynamics 
simulations of a system of Ni = 240 large grains of diame- 
ter Di = 1 cm and density pi — 1.9 g cm -3 , and N s small 
grains with diameter D s and the density p s . We vary D s 
and p s and chose the number N s in such a way that the 
two species occupy a comparable volume, NiDf ~ N S D^. 
The particles are enclosed in a box with a square basis of 
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side length L — 7 cm (see Fig0 with periodic boundary 
conditions in the horizontal directions, so that convec- 
tion is avoided. In order to prevent crystallization some 
particles are randomly glued on the container basis (in 
such a way that no further particle can touch the bottom 
of the container) . 

Two grains interact when in contact via a normal and 
a tangential force. The normal force is given by the so- 
called linear spring-dashpot model, while the tangential 
interaction is implemented by keeping track of the elastic 
shear displacement throughout the lifetime of a contact 
[2^|. The model and the values of its parameters have 
been described in j^, with the value of the viscous co- 
efficient of the normal interaction such that the restitu- 
tion coefficient is e = 0.8. In most of our simulation the 
static friction coefficients are equal for the two species, 
Pu = fJ-is = Pss = 0.4 (pij is the friction coefficient in 
the interaction between a grain of type i and a grain of 
type j, and faj — fiji), but we will also consider the 
case flu ^ p ss to investigate the role of friction in the 
segregation process. 

The system, starting from a random configuration, is 
subject to a tap dynamics up to reach a stationary state. 
Each tap consists of one oscillation of the container ba- 
sis with amplitude A and frequency to, i.e., the bottom 
of the box moves with z(t) — A cos(2iiLot). We checked 
that both A and 10 are important to select the final seg- 
regation state (F is not the only relevant parameter) and 
consider here the case where to — 30 Hz and A is varied. 
A tap is followed by a relaxation time where the system 
comes to rest. A grain is considered to be at rest if its 
kinetic energy becomes smaller than 10 -5 mgd, where 1 
mgd is the energy required to rise it of a distance equal 
to its diameter. All measures are taken when the system 
is at rest and in the stationary part of the tap dynamics. 
Actually, it is known that for small values of T the sys- 
tem dynamics has strong "glassy" features [i^, and 
thus the states attained can be very far from stationarity. 
Here we are away from that region. 

The degree of separation of the binary mixture in the 
stationary state is quantified by the usual vertical segre- 
gation parameter 
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Ah = 2 
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where h p is the average height of particles of species 
p = l,s (h p = l/NpJ^i^i z i> here Zj is the height of 
particle i with respect to the container basis at rest). We 
prepare the system in a random initial state character- 
ized by Ah ~ via a Monte Carlo procedure. When 
subject to a tap dynamics the mixture evolves and the 
segregation parameter changes until a stationary state is 
reached. 

Results- We first describe the dependence of the seg- 
regation parameter, Ah, on the diameter ratio Di/D s 
and on the density ratio pi/p s of the two components 



< 



0.8 
0.4 


-0.4 
-0.8 



• pl/p s =2.0 
pi/p s =1.5 
A pi /ft =1.0 

RBN 

♦ • 

V ♦ 

1.0 1.3 1.6 

Size ratio, Di/D s 



r=2 


, A A £ 




X 


s 


v RBN 






> / 


4 

1 c 


' Di/Ds=1.05 


Di/D s =1.35 
□ D!/D s =1.65 


1 


2 3 4 



Density ratio, pj/p s 



FIG. 2: (color on-line) The left panel shows the segregation 
parameter Ah as a function of the diameter ratio Di / D s of 
the mixture components, for T = 2 and pi/p a = 1, 1.5, 2. By 
increasing Di / D a the system crosses from RBN configurations 
(i.e., Ah > 0) to BN (i.e., Ah < 0). The right panel shows Ah 
as a function of the density ratio pi / p„ for T = 2 and Di/D s = 
1.05, 1.35, 1.65. By increasing pi/p„ the system moves from 
RBN to BN, as the crossover point (where Ah — 0) turns out 
to depend on the size ratio Di/D B . 
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FIG. 3: (color on-line) Ah is plotted as a function of the 
adimensional vibrational acceleration F, in a mixture where 
Di/D„ = 1.25, for the shown values of the density ratio pi/p a . 
As r increases the system moves from a BN to a RBN con- 
figuration. 



for r = 2. Fig. [21 shows Ah as a function of Di/D s for 
different values of the density ratio (left panel) and Ah 
as a function of pi/p s for different values of the size ra- 
tio (right panel). As expected, when the diameter ratio 
grows BN states are favoured with respect to RBN states, 
even though such an effect is mitigated by increasing the 
density ratio pi / p s which shifts the BN to RBN crossover 
to higher values of Di/D s . 

Such a size and density ratio dependence might seem 
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FIG. 4: (color on-line) The "segregation diagram" of the 
mixture in the {pi/ps, Di/D a , F) space. The plot shows the 
regimes where the reverse (RBN) and the usual Brazil nut 
effect (BN) occur. Empty circles are the points where Aft. is 
zero, within 10%, and named "mix" in the caption. The solid 
line separating the areas is given by Eq. J3J. The dashed line 
is the crossover line proposed in Ref.[2j Inset: dependence on 
the adimensional acceleration F of the coefficient a of Eq. (|5J . 



to result from two simple competing effects. The first one 
is a "percolation" effect |4| according to which it is easier 
for the smaller particles to percolate through the voids 
between larger grains and, thus, reach the bottom of the 
container. The percolation effect becomes stronger as the 
size ratio between the components increases, and there- 
fore should describe the size dependence found in Fig. [21 
(left panel). The second effect is buoyancy, according 
to which the system tends to minimize gravitational en- 
ergy, and therefore the species with higher mass density 
is pushed to the bottom of the container. This should 
describe the density dependence of Fig. (right panel) . 

We find, however, that the properties of the external 
forcing have an essential role in selecting the final segre- 
gation state and the overall scenario appears to be richer: 
percolation/buoyancy effects changes with the intensity 
of vibration T (see Fig^J. In Fig. we plot the depen- 
dence of the segregation parameter, Aft, on T for given 
values of Di/D s and pi/p s - unexpectedly, a stronger 
shaking enhances RBN, i.e., as T increases Aft increases 
too. A similar qualitative result was observed also in 
experiments with a continuous shaking dynamics [2l|. 

Fig. ^summarizes these findings in a "segregation di- 
agram" in the (pi/p s , Di/D s , T) space. As expected, the 
BN effect is favored when Di / D s is large and RBN when 
pi/ps grows. The BN to RBN crossover region is depen- 



dent on T: we approximate the BN to RBN crossover 
line Di/D s = f(pi/p s ,T) around pi/p s — 1 with a linear 
function (continuous line in the figure): 



— ~ 1 + cv(r) — 
D s \Ps 
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where the angular coefficient ev(r), shown in the inset of 
FigEl grows monotonically with T. Since a(r) > 0, the 
present results, corresponding to grains with equal fric- 
tion properties (see below for a different case) , point out 
that RBN configurations can be found only if pi/p s > 1: 
i.e., by changing Di/D s there is no way to find RBN 
when pi/p s < 1. In this perspective our simulations 
may explain why in the "original" Brazil Nut Effect ob- 
served during the transportation of nuts of different size 
(but otherwise similar) the larger ones are systematically 
found to rise to the surface. The diagram of Fig. 0] ap- 
pears to be in good agreement with the general features 
of known experiments as those of Ref. [21| (even though 
it is still unclear whether the phenomena of segregation 
under tapping, here considered, and under continuous 
shaking, as in |2lj] , are qualitatively similar) and can help 
clarifying experimental results. 

In Fig. 0] we also plot the BN to RBN crossover line 
found by the theory of Ref.0 (dashed line): such a the- 
ory, approximating the granular mixture under vibration 
as a thermal system, predicts the right qualitative be- 
havior as a function of pi/p s , but it doesn't capture the 
right Di/D s dependence. More elaborated models as- 
suming lack of equipartition between th e sp ecies, as in 
kinetic theories 12] or some simulations 14], or the ex- 
istence of more than one configurational temperature, as 
in a Statistical Mechanics theory of the mixture [lfj, [Fjj , 
may be able to improve on this aspect. 

An other parameter relevant to segregation is fric- 
tion [i"ol [28| which we now consider in order to ex- 
tend the diagram of Fig. 0] We take grains with equal 
sizes and weights, but different friction coefficients; this 
is an interesting situation difficult to be experimentally 
accessed as real grains which differ in frictional prop- 
erties usually also differ in other properties (such as 
mass, Young modulus, etc.). We study a mixture of 
Ni = N-2 = 300 grains of diameter D\ = D-2 = 1 cm 
and density p\ = p2 — 1.9 g cm -3 with friction coeffi- 
cients /in = 0.4, p 2 2 € {0.05,0.01,0.2,0.4,0.6,0.8}, and 
A 4 i2 = niin (/in, /ii2), as in [Fa ]. As the two components 
only differ for their friction, the segregation parameter 
is now defined as Aft = 2(ft,2 — h\)j(h,2 +h%). This mix- 
ture indeed segregates: Fig.0shows that the species with 
higher friction coefficient always rise to the top, as the 
degree of segregation depends on the shaking intensity. 
This can be explained by considering that grains with 
smaller friction can more easily percolate to the bottom 
of the container. 

Conclusions - Our Molecular Dynamics simulations 
are not affected by the presence of air, humidity and (due 
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FIG. 5: (color on-line) The segregation parameter A/i is 
shown as a function of the ratio between the static friction 
coefficient fj,22 / fJ-u of a mixture of grains which differ only for 
their frictional properties. Particles with higher friction coef- 
ficient are always found on the top of the container, althought 
the degree of segregation depends on T. 

to the periodic boundary conditions) convection, and 
should be therefore considered as an ideal, even though 
comprehensive, experiment, well robust to changes in the 
MD model 241. We found that both grain properties, 
such as diameters, densities and friction, and external 
driving properties, such as amplitude and frequency of 
shaking, are important to select the system final segrega- 
tion state. We determined the "segregation diagram" in 
the three parameters space (A/ D Sl pi/p a , T) and derived 
the BN to RBN crossover line A/ A = f(pi/p s ,T). In 
particular, in our model system, a mixture of grains only 
differing in sizes always segregates in a BN configuration, 
explaining why in the original "Brazil nut" problem large 
grains always seat at the top. We also discussed how seg- 
regation is influenced by grains friction by showing, for 
instance, that in a mixture of particles differing only for 
their surface friction by increasing T the smoother grains 
tend to rise to the top, a result easy to be experimentally 
checked. As our results are in agreement with known ex- 
periments and can help their clearer interpretation, our 
comprehensive "segregation diagram" is not explained by 
current theories and necessitates further theoretical and 
experimental investigations. 
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